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(57) Claim 

1-. A method of controlling the rate of oxidation 
of reagents within a mineral heap, comprising: 

taking measurements indicative of physical parameters 
of the heap ; 

analyzing the measurements to predict the fraction of 
reagent remaining in the heap at some future time; 

predicting what fraction of reagent would remain in 
rhe heap at that future time if a measured guantity of one 
-r more control variables were supplied to the heap; and 

then supplying the measured guantity of y'ne of more 
of the control variables which produce the predicted lower 
fraction of reagent at that future time; 

wherein the physical parameters include the 
temperature of the heap and the rate of flow of water 
through the heap; and 

the control variables include heat and water. 

12. A method according to any preceding claim, 
wherein analyzing the measurements involves inputting 
measurements to a computer simulation of the oxygen reaction 
of the reagents within the heap, which uses a mathematical ' 
model of the reaction and takes account of all the measured 
physical parameters. 
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22. A method of designing construction of a heap 
comprising the steps of: 

installing intruments as the heap is constructed to 
take measurements indicative of the physical parameters of 
the heap; 

- installing means as the heap is constructed to supply 
measured quantities of the control variables to the heap; 

determining some of the physical parameters of the 
heap before it is constructed and as it is constructed; 

analyzing the determined values to predict what 
fraction of Reagent would remain in the heap at some future 
time if the value of one or more other physical parameters 
were varied; and 

constructing the heap with said =other physical 
parameters such that the fraction of reagent remaining at 
that future time is reduced or minimized. 
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This invention concerns a method of leaching a heap 
to recover metal, and in particular concerns controlling the 
rate of oxidation of reagents within a mineral heap. in a 
second aspect, the invention concerns a method of designing 
construction oE a heap £or leaching. 

The leaching of metals such as copper, nickel, zinc 
and uranium f rom some ores' depends on the conversion 
(oxidation) of largely insoluble metal sulphides to the much 
more soluble sulphates. The process is usually carried out 
in large heaps of sulphidic material through which water is 
passed. The water dissolves the metal sulphates and passes 
out through the base of the heap where it is 'processed to 
remove the metal of economic interest. The process is 
essentially one of concentration and is generally applied to 
15 low grade ores where conventional hydrometallurgical or 

pyrometallurgical technigues are uneconomic. Crucial to the 
solubilisation process is bacterially catalysed oxidation of 
iron pyrite or other iron containing pyrites. 
* The oxidation of pyrite to ferrous sulphate can be 

.20 described by the eguation. 



10 



? «. * 



FeS 2 + jo 2 *■ H 2 0 " » FeS0 4 + H 2 S0 4 O) 

The oxidation of ferrous sulphate to ferric sulphate 
.25 can in turn be described by, 

FeS0 4 + ^0 2 + ^H 2 S0 4 > iFe 2 (S0 4 ) 3 + ^0 (b) 



30 



Under the correct chemical conditions ferric sulphate 
will oxidise more pyrite according to the equation 

FeS 2 + 7FE(S0 4 ) 3 + 8H 2 0 > 15 FES0 4 + 8H 2 S0 4 ( c ) 

' r ^ X he last two equations describe an autocatyli t ic 

process which is also catalysed by. the intervention of 
microorganisms which speed up the second reaction. 
35 Ferric sulphate will also convert other metallic 

sulphides to sulphates according to the general equation, 
MS + 2FE 2 (S0 4 ) 3 + 2H 2 0 + 0 2 > MS0 4 + 4FeS0 4 + 2H 2 Sp 4 (dj 
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It is also known that the extraction efficiency of 
gold from many auriferous sulphide ores is greatly increased 
if the ore is first oxidised. .This can be effected ^in a 
heap, where the oxidation process can be described by 

5 equations (b) and (c) ' above ^ 

In an industrial ^heap, which may measure tens to 
hundreds of meters in height, the oxidation rate is limited 
by the rate at which oxygen can be supplied to the sites of 
unreacted pyrite. 

10 Initially oxygen will enter the heap from the 

atmosphere/heap interface by the process of molecular 
diffusion. At this time, the oxidation reaction is confined 
to the surface layers of the heap. The heat released from 
the pyritic reaction will-induce a convective air current 

15- within the heap "and in so doing will increase the rate of 
atmospheric oxygen transport into the heap and permit" 
.oxidation throughout the heap. The oxygen supply to the 
heap is believed to be the limiting factor of the oxidation 
reaction. 

20 According to a first aspect of the present invention, 

there is provided a method of controlling the rate of 
oxidation of reagents within a mineral heap, comprising: 

taking measurements indicative of physical parameters 
of the heap; 

^25 analyzing the measurements to predict the fraction of 

reacent remaining in the heap at some future time; 
* predicting what fraction of reagent would remain in 

the heap at that future time if a measured quantity of one 
or more control variables were supplied to the heap; and 
30 then supplying the measured quantity of one of more 

of the control variables which produce the predicted lower 
fraction of reagent at that future time; 

wherein the physical parameters include the 
temperature of the heap and the rate of flow of water 
3 5 through the heap; and 

the control variables include heat and water. 
Taking measurements indicative of the temperature of 
the heap includes taking the temperature at each point of an 
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array distributed throughout the heap or part of the heap. 
It also includes non-invasive techniques of obtaining 
measurements indicative of temperature, such as those 
employing infra-red detection. 

5 Taking measurements indicative of the rate of flow of 

water through the heap (irrigation rate) includes measuring 
the net volume of water flowing from the bottom of the heap, 
and measuring the rates of flow in particular directions at 
points throughout the heap or parts of the heap. 

10 Measurements may also be taken indicative of oxygen 

concentration, pyrihe concentration, air pressure, particle 
size and the geometry of the heap. All these variables have 
an effect on the convective and diffusive transport of 
oxygen into the heap. Measurements may also be taken to 

15 determine gas permeability and diffusion coefficients in the 
heap. 

Analyzing the measurements includes weighting, 

[] c averaging or otherwise adjusting the measurements, and 
preferably includes inputting the measurements to a 

t 20 simulation of the oxidation reaction of the reagents within 
the heap. The simulation would usually be performed on a 
computer using a mathematical model of the oxidation 
reaction. The simulation will advantageously take into 
account all the measured physical parameters of the heap. 

25 Supplying a measured quantity of one or more of the 

control variables changes the amount of heat in the heap. 
The heat in the heap may be increased by direct means, such 
as electric heating elements positioned in an array 
throughout the heap, or a selected portion of the heap such 

30 as the bottom 1 metre; it may also be supplied by 

irradiating the heap, say, with microwave radiation. 

l' M Heat may be indirectly supplied to the heap by 

increasing the convection rate of air or oxygen in the heap, 
for instance by injecting air or oxygen through pipes into 

35 the interior of the heap. Air or oxygen may, of course/ be 
included in the list of control variables. 

Removing heat from the heap can be achieved by 
increasing the flow of water through the heap, such as by 
pouring or sprinkling water onto the heap, or injecting 
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water into the heap from an array of nozzles arranged ^ 
throughout the heap or a portion of it. 

Water flowing through the heap cools it in two ways. 
Firstly, it removes heat directly by absorbing thermal heat 
and latent heat. Secondly, it removes heat by obstructing 
the interparticular space in the heap and slowing convection. 

It is important that both heat supply and heat 
removal are controlled since too much heating could raise 
the temperature within the heap to a point where the micro- 
organisms cease to be effective as catalysts, but too little 
heat will slow down the oxidation reaction. 

Advantageously, local environmental and economic 
factors are also taken into account; environmental factors 
include the precipitation rate, and the evaporation rate. 
15 Ideally, measurements indicative of these factors will be 
taken and supplied to the simulation. 

The analyzing and predicting steps are preferably 
repeated at regular intervals, under the control of an 
operator. 

e 9 0 Optionally, the method includes the step of 

installing intruments in the heap to take measurements 

V* indicative of the physical parameters. Also, a step of 

installing means to effect the supply of the control 
variables may be included. 
25 According to a second aspect of the present invention 

as currently envisaged, there is provided a method of 
designing construction of a heap comprising the steps of: 

installing intruments as the heap is constructed to 
take measurements indicative of the physical parameters of 

30 the heap; 

installing means as the heap is constructed to supply 
measured quantities of the control variables to the heap; 

determining some of the physical parameters of the 
heap before it is constructed and as it is constructed; 
35 analyzing the determined values to predict what 

fraction of reagent would remain in the heap at some future 
time if the value of one or more other physical parameters 
were varied; and 
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constructing the heap with said other physical 
parameters such that the fraction of reagent remaining at 
that future time is reduced or minimized. 

The invention will now be described by way -Of example 
only with reference to the accompanying drawing, which is a 
graph showing the behaviour of a heap under various 

circumstances. 

instruments are applied to a heap to measure the gas 
concentration, pyrite concentration, temperature, irrigation 
rate (flow of water through the heap), air pressure, and 
particle size; the geometry of the heap is also measured. 

Measuring these variables serves two purposes. 
Firstly, it monitors changes in the heap as the oxidation 
proceeds; and secondly, it allows values to be generated 
which are fed to a computer model to obtain predictions of 
the future behaviour of heap. An operator will decide at 
what intervals it is appropriate to feed parameters to the 

computer model . 

The computer model is able to simulate the oxidation 
reaction in the heap based on the information received about 
the physical parameters. The prediction is usually 
expressed in terms indicative of the fraction of reagent 
left in the heap when compared to some reference date. The 
fraction of reagents left within the heap may be indirectly 
expressed in terms of some convenient statistic which, for 
the purpose of convenience, is called the global oxidation 
rate or GOR. 

The reference date may be selected in any convenient 
manner. For instance, the reference date could be the date 
metal recovery commenced, or the date of the analyzing 
step. The date may be changed every time the analyzing step 
is performed, or may remain static in time throughout all or 
part of the metal recovery process. 

In general, the operator will decide from the 
information available to him (about the physical parameters, 
the environmental factors and any relevant economic 
constraints), whether it is appropriate to increase the 
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temperature within the heap or to reduce the temperature 
within the heap. If the operator decides an increase in 
temperature is called for, he will estimate what quantity of 
heat to supply to the heap, whether to supply it directly, 
or indirectly by reducing irrigation or increasing air or 
oxygen convection, and how it is to be distributed within 
the heap. He will then model the behaviour of the heap 
using a computer model to predict the likely effect of the 
change estimated. When an estimated change produces an 
acceptable prediction, the operator will cause that change 
to be effected. 

If the operator decides a reduction in the temperture 
in the heap is called for, he may do this by reducing the 
amount of heat supplied, by supplying water to the heap, or 
by reducing the supply of air or oxygen to the heap. He 
will then go through the estimation and modelling procedure 
until a desired result is predicted, and then effect the 
appropriate change. 

The computer model is based on a deterministic 
formulation of multiphase flow through porous media as 
follows: 

The mass balance equations for each phase 

a = g, J, s 

(g = gas(air), I = liquid(water ) , s = solid (heap ore}) 
are given by 

V.(P a v a ) = w a (1) 

where is the bulk density of the a phase and is related 
to the intrinsic density 0 by 



where is the volume fraction of the a phase such that 
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30 



The quantity W 0 on the right hand side of (1) is an external 

source/sink of the a phase. For the liquid and gas phases 
Darcy's Law is assumed: 

e v« = - ^IS (v P *'p°gcJ ■ (4) 

where K is the saturated permeability; k ra is the relative 
permeability of the a phase which is some function of e a ; p^ 
is the pressure of the a phase; U a is the viscosity of the a 
phase; g the acceleration due to gravity and e z is the 
unit vector pointing verticaly upwards. We assume that 

V s = 0 (5) 

and that P and E g are constant despite the loss of 

mass to the 5 solid phase due to the oxidation of the pyrite. 

The mass balance equation for a species iG{l., . . . ,n Q } 
in the a phase is given by 



3tJ i a . . ..a a . t T a = , T . * + G a (6) 

p a IF 



;?0 -where u. a is the mass fraction of the species i in the a 
phase; J- a is the nonadvective flux vector for species i 
in the a phase, I. a is the diffusive mass fluxes across 
the phase interfaces and g/* is an external source/sink 
term of the species i in the a phase. 
25 in the gas phase, we consider only the oxygen specie 

so that 



where w 9 is the mass fraction of oxygen in the gas phase 
and the subscript i has been dropped since only one specie 
is considered. G 9 is the injected oxygen. The 
nonadvective flux vector J g is assumed to be described by 
the Fickian law 

35 j9 = _ p D v„9 (8) 

u g g 
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where D is the coefficient of diffusion of oxygen in the 
porous medium. It should be noted that we have ignored 
dispersion which could be significant when convection gets 
going. This exclusion should not be seen as a restriction 
5 to the model as it is straight forward to insert it at a 

later stage. 

The pyrite is fixed in the solid phase so that {6) 
reduces to 

do: 5 „ T s (9) 

where * s is the mass fraction of pyrite in the solid phase 
and again the subscript i has been dropped since there is 
only one specie considered in the solid phase. Note that 
15 here we exclude any external source/sink terms. 

The energy eguation for the porous medim is given by 



* * * 

3. 3 



, ,?0 



25 



30 



35 



where T is the temperature; c a is the specific heat capacity 
of the <* phase; D h is the thermal conductivity in the 
porous medium (all phases), S is the heat produced by the 
oxidation reaction and S exfc is the externally applied heat. 

Both p s and P 1 are assumed to be constant whereas 
the intrinsic gas density is assumed to be a function of 
temperature through 

0 9 = p g (1 - aT) 

0 

where p% is the atmospheric air density and a is the 
coefficient of thermal expansion. 

The boundary conditions are obtained by setting the 
oxygen concentration, temperature and air pressure at their 
atmospheric values at the atmosphere/heap interface. On the 
upper boundary of the heap is imposed a flux condition. 

_ (Vp +P 'ge ).n = -o'gCQg+Qj) (12 > 
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which induces vertical water infiltration which is assumed 
to be drained at the base of the heap. In (12) Q R is the 
net input of water by rainfall/evaporation and the 
input by artificial means, 

5 Upon incorporation of the constitutive relations of 

(2), (3), (4) and (11) the 2 equations of (I) a = g,/; 
equations (7) and (9) and the heat equation (10) constitute 
a system of 5 equations for the 5 variables p g , Pj, 
co g , w s and T. We define the "natural" operation of the 

10 above heap system when there exists no artifical 
inputs/outputs Qj, S fc# G g , W a ( *=f,g). Under 
these natural conditions, diffusion will initially be the 
only means of transport of oxygen to reaction sites within 
the heap. As the heat of oxidation builds up convection 

15 will develop and raise the GOR. The extent to which both 
diffusion and convection are effective in transporting 
oxygen into the heap will depend on the physical parameters 
1 *« which define the particular heap system under 

consideration. The physical parameters which are particular 

20 to each heap are contained in a set 
;f P H{0 s ,K / e s ,a ) S i £=0 ,D g ,c s ,D h ,k ra (£ a )(a=I < g),Q R , 

4 -t 

L _ 

i g (^ g ^ S ,T),I S (a 3 g ^ S ,T) ,S(a J g .^ S ,T)L 1 L m . 

The interdependent quantities I g , I s , S define the 

25 oxidation rates and heat rate of reactions and are defined 

by a set of parameters which must also be included in the 

set P. These reaction parameters will depend largely on the 

microscopic structure of the medium. Also in the set P we 

*** include a set of parameters L, , . . . , L which define 

l m 

30 the geometry of the heap. For the example of a cylindrical 
heap this set will simply contain the radius R Q and the 
- height Z Q . 

In a preferred form of the second aspect of the 
invention, the heap is constructed with instrumentation for 
35 monitoring. Simulations using the computer model are 

performed prior to construction of the heap in order to 
decide on an optional set of heap parameters defined in P. 
This may entail modificaion of run of mine material. Then 
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data ir>»-« are fchen fed as initial 

data into the computer model which simulates the 

7 1Utl0narY behaviour of the operation and after a series 
--- the operator decides on a course of ac on w 

ool L I™ G ° R - S ° me dSCisi - facilities 

could be incorporated in the model. The course of t - 
effected , y the special variation Qf _ » 

" -put entities Qj . s ^ ^ 

T«e computer model is based on a finite differ 

th : ~— *™" e b r r 

set of it , ° P ° E Speci »^ Prefers. Given . 

set of mitral data fer the variables p , „ g 

» - and T (which ace obtained by insitu measurement) ^he 

:::::: run to sim — ■ ■>«*«- 

interval. At the end of this time interval tho 
Progr,™ computes the current GOR If "tar. 

in the GOR the operator win recun " ^ ln """ 

<=ha„ ges i„ one or more of the artilicra! " '"" lel 

S w r , , artificial paramaters 0 

nSt'til = SUCh that G ° R iS --eased in the' 

next: time lntprval r.p *-v,~ , c 
j-uttiLvai ot the heap oDershinn 0 

features could he incorporated in the lael ^-f '""^ 
rs repeated sequentially until either The ™ """^^ 

. is depicted or the heap ^tZT^T" ^ >"> 
Although the invention has been ITcrZ^T"" ' 
reference to preferred embodiments, it should he 
^at it could be performed in many other 'V"™""- 
the mathematical desctiot ln „ « t I instance, 
include the chemica s !c Is „ £ « also 

The transport of a „ y ehemica! P e C U oTl "1"' ^ 
tof. i = l „ specie, of mass fraction 

, 6) wi / h 15 9 ove ™ed by n , equations of 

7 tor a -f to account for 
homogeneous reactions of fh Q - 

Here « b^^ IT "o^ 
artrfrcral i nputs and the v „ iablM „, «' 
are added to the set- of ™ - 1 • • • • . n. 

set of monitored variables of the heap. 
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Referring now to Figure 1: 

(1) shows the behaviour of the heap with no 
artificial heating and no irrigation. 

(2) shows the behaviour of the heap with continuous 
artificial heating to the bottom 1 metre of the heap at 
100 W/m 3 and no irrigation. This graph shows an initial 
improvement of the GOR over (1) but overheating eventually 
retards the GOR. 

(3) shows the behaviour of the heap with an 
irrigation rate of 0.1 m per hour and no heating. In this 
case, heat has been removed from the heap and convection has 
been reduced which results in a decreased GOR, 

(4) shows the behaviour of the heap with an 
irrigation rate of 0.01 m/hour and heat supplied to the 
bottom 1 metre of the heap at a rate of 100 W/m . In this 
case, significant air flow has been induced which results in 
overheating and the reaction switching off. However, it is 
worth noting that there is some improvement over the 
situation where only the heat was applied (2) and the 
situation where only the irrigation was supplied (3). 

(5) shows the behaviour of the heat with an 

—3 

irrigation rate of 10 m/hour with no artificial heating, 
and this shows a considerable improvement over the natural 
state. 

(6) shows the results obtained from a process 
embodying the present invention, where the irrigation rate 
is tailored at regular time intervals. This shows 
considerable improvements over all previous results. 

(7) also shows the result of a process embodying the 
present invention. Here, the bottom 1 metre of the heap was 
artificially heated according to a predicted amount and a 
constant rate of artificial irrigation of 0.1m per hour was 
applied. 
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THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS: 

1. A method of controlling the rate of oxidation 
of reagents within e mineral heap, comprising: 

taking measurements indicative of physical parameters 
5 of the heap; 

analyzing the measurements to predict the fraction of 
reagent remaining in the heap at some future time; 

predicting what fraction of reagent would remain in 
the heap at that future time if a measured quantity of one 
10 or more control variables were supplied to the heap; and 

then supplying the measured quantity of one of more 
of the control variables which produce the predicted lower 
fraction ur reagent at; that future time; 

wherein the physical parameters include the 
15 temperature of the heap and the rate of flow of water 
through the heap; and 

the control variables include heat and water. 

2. A method according to claim 1, wherein taking 
measurements indicative of the temperature of the heap 

20 involves taking the temperature at each point of an array 
distributed thoughout the heap. 

3. A method according to claim 1, wherein taking 
measurements indicative of the temperature of the heap 
involves infra red detection. 

25 4. A method according to any preceding claim, 

wherein taking measurements indicative of the irrigation 
rate involves measuring the net volume of water flowing from 
the jottpra of the heap. 

r: , : -,- :-^ r:r:: ' : ""' 5. A method according to any one of claims 1 to 3, 

30 wherein taking measurements indicative of the irrigation 
rate involves measuring the races of flow particular 
directions of points throughout the heap. 

6. A method according to any preceding claim 
wherein measurements indicative of oxygen concentration are 

35 taken. 
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7. A method according to any preceding claim, 
wherein measurements indicative of pyrite concentration are 
taken. 

8. A method according to any preceding claim, 

5 wherein measurements indicative of air pressure are taken. 

9. A method according to any preceding claim, 
wherein measurements indicative of particle size are taken. 

10 . A method according to any preceding claim, 
wherein measurements are taken of the geometry of the heap. 

10 11. A method according to any preceding claim, 

wherein measurements are taken to determine gas permeability 
and diffusion coefficients in the heap. 

12. A method according to any preceding claim, 
wherein analyzing the measurements involves inputting 

15 measurements to a computer simulation of the oxygen reaction 
of the reagents within the heap, which uses a mathematical 
model of the reaction and takes account of all the measured 
physical parameters. 

13. A method according to any preceding claim, 

20 wherein the heat in the heap is increased by direct means. 

14. A method according to any one of claims 1 to 
12, wherein heat is indirectly supplied to the heap by 
increasing the convection rate of air or oxygen in the heap. 

15. A method according to any preceding claim, 
25 wherein the control variables include air. 

16. A method according to any preceding claim, 
wherein the control variables include oxygen. 

17. A method according to any preceding claim, 
wherein removing heat from the heap is achieved by 

30 increasing the flow of water through the heap. 

18. A method according to any preceding claim, 
wherein the analyzing and predicting steps are repeated at 
regular intervals under control of an operator. 

19 . A method according to any preceding claim, 
35 wherein the method includes the step of installing 

instruments in the heap to take measurements indicative of 
the physical parameters. 
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-20. A method according to any preceding claim, 
wherein the method includes the step of installing means to 
effect the supply of the control variables. 

21. A method of controlling the rate of oxidation 
of reagents within a mineral heap substantially as 
hereinbefore described with reference to the accompanying 
drawings. 

22. A method of designing construction of a heap 
comprising the steps of: 

installing intruments as the heap is constructed to 
take measurements indicative of the physical parameters of 
the heap; 

installing means as the heap is constructed to supply 
measured quantities of the control variables to the heap; 

determining some of the physical parameters of the 
heap before it is constructed and as it is constructed; 

analyzing the determined values to predict what 
fraction of reagent would remain in the heap at some future 
time if the value of one or more other physical parameters 
20 were varied; and 

constructing the heap with said other physical 
parameters such that the fraction of reagent remaining at 
that future time is reduced or minimized. 
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DATED 9th day of August, 19 9 0 
AUSTRALIAN NUCLEAR SCIENCE AND TF CHNPLOCV f>T>GAWTRATTOH 

By their Patent Attorneys 
GRIFFITH HACK & CO. 
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